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Abstract 





In this presentation, a discontinuous Galerkin augmented electric field 
integral equation (DG-AEFIE) method based on the domain 
decomposition is proposed for full-wave solution of multiscale targets. 
The common surface integral equation based discontinuous Galerkin 
method allowing both conformal and nonconformal discretizations 
for multiscale structures suffers from the low-frequency breakdown. 
By augmenting DG-EFIE with current continuity equation, the 
proposed scheme can alleviate the low-frequency breakdown. In the 
augmented system, the EFIE and the current continuity equation are 
discretizated by using hybrid basis functions including RWG and half 
RWG basis functions. Since the half RWG basis is not divergence- 
conforming, line charge degrees of freedom on the adjoining edge are 
introduced in this work. It is observed that the resulting linear system 
is well-conditioned at low frequencies, which leads to a rapid 
convergence over wide frequency band. 


Keywords: Discontinuous Galerkin, electric field integral equation, 
low-frequency breakdown, multiscale problems. 
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€ Domain Decomposition : Divide and Conquer 
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€ Non-Conformal Mesh 
€ RWG basis functions 
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Background-(DG-EFIE) 


@ Impedance Matrix of EFIE (Half RWG) 
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€ Non-conformal discretization 
© Different basis functions can be seamlessly integrated 
9 lp-adaption: higher order bases (p refinement), 


refinement of the mesh (h refinement) 
€ Low frequency breakdown 
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@ Augmented electric filed integral equation (AEFIE) 
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€ DG-AEFIE 
© DD of PEC surface 


9 Choice of basis functions 
Half RWG basis functions + RWG basis functions 
Surface charge basis functions + Line charge basis functions 
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DG-AEFIE Formulation 








€ Surface current expansion 
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DG-AEFIE Formulation 
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9 DG-AEFIE with charge neutrality 








9 Linear system 
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@ DG-AEFIE with preconditioner 


9 Left preconditioned linear system 
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© Approximation of impedance matrix 
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Sphere-Accuracy 





€ Plane wave with amplitude of 1V/m , along the -z axis, 


© Generalized conjugate residual (G 


CR) and LU direct method, 


© Relative residual error (RSS) of 10-8. 
Discretization of Sphere: 1458 triangles, 731 nodes and 2187 edges. 
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Condition number of the impedance matrix and 
relative error of current as the stabilization parameter 
B varies at 300 MHz. 
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Sphere-Conditioning 








Discretization of Sphere: 1458 triangles, 731 nodes and 2187 edges. 
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of five linear systems. 





Surface current density calculated by using EFIE-LS 
and DG-AEFIE (electrically size 7.1x10 ^A) 
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Cone-Nonconformal Mesh 
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Chip-Nonconformal Mesh 





€ Plane wave excitation 
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€ Voltage source excitation 
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Conclusion 





€ DG-AEFIE can handle complicate structures with 
nonconformal mesh, 


€ DG-AEFIE is well-conditioned from the low 
frequencies to high frequencies, 


€ DG-AEFIE allows local refinement of the mesh. 
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